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We show that the relativistic dissipative hydrodynamic equation derived from the relativistic 
Boltzmann equation by the renormalization-group method uniquely leads to the one in the energy 
frame proposed by Landau and Lifshitz, provided that the macroscopic-frame vector, which defines 
the local rest frame of the fluid velocity, is independent of the momenta of constituent particles, as 
it should. We argue that the relativistic hydrodynamic equations for viscous fluids must be defined 
on the energy frame if it is consistent with the underlying relativistic kinetic equation. 
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I. INTRODUCTION 

Theory of relativistic hydrodynamics for viscous flu- 
ids is a powerful mean for analyzing the long-wavelength 
and low-frequency dynamics of many body systems, such 
as the hot and/or dense quark-gluon or hadronic matter 
created by relativistic heavy-ion collisions [1-4] and also 
various high-energy astrophysical phenomena [5, 6]. 

However, the relativistic dissipative hydrodynamic 
equation is still under debate in the fundamental level. 
Indeed, the following three problems may be noted: (A) 
There are ambiguities in the definition of the fluid veloc- 
ity [7, 8], (B) In the Eckart (particle) frame, there arises 
an unphysical instabilities of the equilibrium state [9], 
and (C) the so-called first-order equations lack in causal- 
ity [10-12]. 

As is easily thought of, it is a legitimate and natural 
way for obtaining the proper relativistic hydrodynamic 
equation, to start with the relativistic Boltzmann equa- 
tion (RBE) which is Lorentz invariant and does not have 
stability nor causality problems [10]. The problem is 
how to obtain an asymptotic dynamics in the far-infrared 
long-wavelength limit of the RBE, and several reduction 
methods [13-15] have been developed with some but not 
a compete success. 

Recently, for examining the first two problems, (A) and 
(B), Ohnishi and the present authors [16, 17] applied the 
renormalization-group (RG) method [18-21] as a power- 
ful reduction theory of dynamics to the RBE and suc- 
ceeded in deriving a generic form of relativistic dissipa- 
tive hydrodynamic equations with the stable equilibrium 
state. A key ingredient in the derivation was introduction 
of a time-like Lorentz-covariant vector a M , with a > 0: 
a M specifies the macroscopic and covariant coordinate 
system where the local rest frame of the fluid velocity 
it M is defined, and is called the macroscopic-frame vec- 
tor. a M could depend on the momenta p^ of constituent 
particles of the system as well as the space-time coordi- 
nate x^. In fact, Ohnishi and the present authors [16, 17] 
could manage to derive the hydrodynamic equation in 
the Eckart (particle) frame only by making have a 



dependence. In retrospect, the possible momentum de- 
pendence of a M , however, may not be legitimate for a M to 
play a macroscopic-frame vector, because it means that 
the macroscopic space-time is defined for respective par- 
ticle states with a definite energy-momentum, and may 
lead to a difficulty in the physical interpretation of the 
space and time in which the hydrodynamics is defined. 
Thus, we are lead to require that a M should be indepen- 
dent of p* 1 and time-like vector with the Lorentz covari- 
ance. 

In this paper, taking this requirement from the out- 
set, we shall examine the outcomes and reach a signifi- 
cant conclusion that the relativistic hydrodynamics con- 
sistent with the underlying kinetic equation should be 
uniquely the one in the Landau-Lifshitz (energy) frame, 
but not in other frames. This means also that the so- 
called matching conditions for selecting the energy frame 
but not other frames [10] is uniquely derived from the 
underlying kinetic equation [22]. 

This paper is organized as follows: After a brief ac- 
count of the basic properties of the relativistic Boltzmann 
equation (RBE), we show that that is independent of 
the momentum p M must be naturally proportional to the 
fluid velocity u M ; a M = bu^. Furthermore, we show that 
the "normalization" factor b can be made unity with- 
out loss of generality. Then, we apply the RG method 
to derive the relativistic hydrodynamics from the RBE 



with = 



and show that the resulting hydrody- 



namic equation is uniquely in the energy frame. We also 
clarify the meaning of the matching condition in terms 
of the inner product for the distribution functions. The 
last section is devoted to a summary and concluding re- 
marks on the cases of multi-component systems and the 
extended thermodynamics. 



II. RELATIVISTIC BOLTZMANN EQUATION 

To make the presentation self-contained, we first sum- 
marize the basic facts about the relativistic Boltzmann 
equation (RBE) [10] very briefly. The RBE is an evo- 
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lution equation of the one-particle distribution function 
f p (x) denned in the phase space (x, p), 

P»dJ p (x) = C[f] p (x), (III) 

with p^ being the four-momentum of the on-shell parti- 
cle, i.e., p^p p = p 2 = m 2 and p° > 0. The C[f] p (x) in 
the right hand side denotes the collision operator, 

C[f] P (x) = ^ i i Yl 4^(P,Pib2,p 3 ) 

PI 1 P2 2 P3 3 

x (fp2 ( x ).fp3 (x) - .f P (x).f Pl (x)j , (II.2) 

with ui(p , pi \p2 , Ps) being the transition probability, 
which has the symmetry property uj(p,p\\p 2 ,p 3 ) — 

W(P2,P3|P,P1) = W(pi,p|p 3 ,p 2 ) = Uj(p 3 ,p 2 \pi,p) 

and respects the energy-momentum conservation 

u(P,Pl\P2,P3) <* S 4 (p + p 1 -P2 -Pz). 

Thanks to the above symmetry property and the 
energy- momentum conservation, the function <p p (x) = 
a(x) +p^ b p (x) is found to be a collision invariant; 

Ei^)W = o, (ii.3) 
p 

where a{x) and b^(x) being arbitrary functions of x. On 
account of Eq. (II. 3), we have formally the balance 
equations for the energy-momentum tensor T^ v (x) = 
J2 P "tt P^ P v fp( x ) an d the particle current N^(x) = 
E P % P 11 f P {x) as follows; d v T^(x) = and d^(x) = 
0, respectively. Despite the appearance, these equations 
do not have any dynamical information before the evolu- 
tion of f p (x) has been obtained from Eq. (II. 1). 

The entropy current is defined by S^(x) = 
~E P pT>P^fp( x )(M 27T ) 3 f p (x) - 1), which is generically 
non-conserving; 

d»S»{x) = -^i« ln(2 7r) 3 /p(^), (H.4) 
p 

due to Eq.(II.l). One sees that becomes a conserved 
quantity only when ln(2 7r) 3 f p (x) is a collision invariant, 
ln(2 7r) 3 f p (x) = a(x)+p^ b p (x). Accordingly, an entropy- 
conserving distribution function can be cast into the form 
of the local equilibrium distribution function (the Jiittncr 
function [23]) in terms of the local temperature T(x), 
chemical potential fi(x), and time-like fluid velocity 1^(21) 
as f p (x) = (2^)- 3 exp[( M (x) - p^(x))/T(x)] = f^(x), 
with u^(x)u p (x) = 1 and u°(x) > 0. We also note that 
the collision operator identically vanishes for the local 
equilibrium distribution function f p q (x); 

C[r\(x)=0, (11.5) 

due to the energy-momentum conservation in the colli- 
sion process. 



III. MACROSCOPIC-FRAME VECTOR 

To implement the task to solve the RBE in the hydro- 
dynamic regime, it is customary [10, 16, 17] to introduce 
a time- like Lorentz vector a M , with a > 0. We call 
the macroscopic-frame vector, which defines the covari- 
ant and macroscopic coordinate system (t, a^) from the 
space-time coordinate x^ as 

dr = a^dx p , (III.l) 
do* = {g»» -^f-)dx v . (III.2) 

It is noted that (r, a^) is the so-called local Lorentz frame 
when a* is dependent on x*. 

We now show that a* must be proportional to the fluid 
velocity provided that a* should be independent of 
the momentum p^ and time-like vector with the Lorentz 
covariance. Noting that u* and d* are the only available 
Lorentz vectors at hand, we see that the generic form of 
a Lorentz-covariant vector reads 

a* = A 1 u" + A 2 d»T + A 3 d^fj, + A A u v d u u^, (111.3) 

where Ai with i = 1, 2, 3, 4 is an arbitrary Lorentz-scalar 
function of the temperature and the chemical potential; 
Ai = Ai(T, /i). Now the derivative <9 M can be decom- 
posed into the time-like and space-like components as 
gn = u v u v q v + ( g »v _ u n u v)d v = u^D + V*, where 
D = u u d v and V = A» v d v with A* v = (g» v - u* u v ) 
being the projection operator onto a space- like vector or- 
thogonal to the time- like vector u* . Then Eq.(III.3) is 
also decomposed as 

a* = (Ai + A 2 DT + A 3 Dp) u* 

+ A 2 WT + A 3 V f "p + A 4 Du f ". (III.4) 

But since a* should be time-like, we have 

= (Ai +A 2 DT + A 3 Dp) u» = b(T, /z) u*, (III.5) 

with b(T, /i) > 0. Here, we have used the fact that V^T, 
V M /i, and Du* are space-like because A* u V„ = V M and 
A^ Du v = Du^. 

With the use of this macroscopic-frame vector, the co- 
variant and macroscopic coordinate system (t, a*) reads 
dr = b(T, p)u^dx tl and da* = A* v dx v . In fact, the 
"normalization" factor b(T, n) is a redundant degree of 
freedom for the dynamics because the b(T, /i) can be al- 
ways made unity by converting r into the new temporal 
coordinate r' as dr' = b(T, dr = dx^. From now 
on, we thus set 

a» = ur (III.6) 

In terms of the new coordinates (r, er^), Eq.(II.l) is 
rewritten as 

^.fp(r, a) = ^C[f] p (r, a) 

_ £ _Lp. V / p (r, a), (III.7) 
n ■ ii. 
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where d/dr = D and d/dcr^ = V M . Here the small pa- 
rameter s, which will be set back to unity eventually, 
represents the non-uniformity of space. This seemingly 
mere rewrite of the equation reflects a physical assump- 
tion that only the spatial inhomogeneity is the origin of 
the dissipation. We note that the RG method applied to 
the non-relativistic case with this assumption successfully 
leads to the Navier-Stokes equation [25], which means 
that the present approach [16, 17] is a covariantization of 
the non-relativistic case. 



IV. REDUCTION TO HYDRODYNAMIC 
EQUATION 

In this section, we derive the relativistic hydrodynam- 
ics through solving the converted Boltzmann equation 
(ni.7) in the RG method. 

We should note that some of formulas presented in this 
section may be found in the previous papers [16, 17] in 
a different context and hence in a different order. Thus, 
the presentation will be a brief reorganization of them; 
we refer to Ref. [17] for the detailed reasoning. 



A. Hydrodynamics from relativistic Boltzmann 
equation by RG method 

In the RG method developed in Ref.'s [16, 17, 19, 20], 
we first try to obtain the perturbative solution f p to 
Eq. (III. 7) around the arbitrary initial time t = tq with 
the initial value / p (t , a); f p (r = t , cr; r ) = / p (t , a). 
Note that the solution depends on the initial time r at 
which f p (r — tq, a; To) is supposed to be on an exact so- 
lution. We expand the initial value as well as the solution 
with respect to e as follows: 

oo 

f p (r, a- T ) = £l °> T ")' ( IV -!) 

i=0 
oo 

/ p (r ,<7) = JV/^fo.")- (IV-2) 

i=0 

The zeroth-order equation is 

^/i 0) (r, a; r ) = j-^ C\J<-%(t, a; r ). (IV.3) 

Since we are interested in the slow motion realized 
asymptotically for r — > oo, we take the stationary so- 
lution satisfying dfp°\r, a; t )/8t = 0, which demands 
that 



C[/(°)] p (t, a; r ) = 0, Ma. 



(IV.4) 



Such a solution is given by a local equilibrium distribu- 
tion function, i.e., the Jiittner distribution function, 



/i 0) (r, a- to) 



1 



exp 



(2tt) 3 

= to), 



fj,(a; to) -jfu^a; t ) 



T{a; to) 



(IV.5) 



with u^{<j] To)u^(a\ To) = 1. We note that the inte- 
gration constants T(a; To), /Lt(<r; To), and u fJ _(a; To) arc 
independent of r but may depend on t as well as a. 
The first-order equation is given by 

y^(r)=Y,A pg f^(r) + F p , (IV.6) 



dr 



where 



p-u p 



(IV.7) 



with A pq being a matrix element of the linearized collision 
operator A, 



{A) pq — Apq 



P-U df q 



(IV.8) 



We now show that the linearized collision operator A 
has remarkable spectral properties [16, 17], which are es- 
sential for the following analysis. We first convert A to 
another linear operator, 



(IV.9) 



with (f eq ) pq = fp q S pq , which is diagonal. Next, an inner 
product is introduced for arbitrary nonzero vectors ip and 
ip as 



which satisfies the positivc-dcfinitcncss of the norm as 
(<P, f) = E P pV (P ■ u) /° q {ifipf > 0, for ip p jt 0, due to 
(p ■ u) > because of both p M and u M being time-like 
vectors. 

Then, as is shown in Ref.'s [16, 17], the linearized col- 
lision operator is semi-negative definite and has five zero 
modes given by 



Vo P = 



p^ for a = fi - 
1 x m for a = 4. 







(IV.ll) 



The functional subspace spanned by the five zero modes 
is called the Pq space and the projection operator to it 



is denoted 



by P ; [Poip] p = VopVoap (<Po > V0> 



where 



T)oa0 ls ^ ne mverse matrix of the the Po-space metric ma- 
trix ^ defined by rfc = (ipg , tp$). We also call the 
complement to the Po space the Qo space and introduce 
Qo = 1 — Pq. In the following, we also use the modified 
projection operators defined by P = f cq -Po(/ cq ) _1 an d 
Oo = / eq Qo(/° q )- 1 - 

Then, the perturbative solution up to the second order 

reads / p (t, a; t ) = jf } + ef™ + s 2 /^ + 0(e 3 ), where 



= (t - t )P F - A- X Q F, 



(IV.12) 
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and we refer to Ref. [17] for a lengthy formula of fp , for 
the sake of space. We note that secular terms are present 
in jp^ (i = 1, 2, • ■ •) , which is caused by the zero modes 
of the linearized collision operator A. 

Remarks are in order here: In the usual approaches 
[10, 17] to derive hydrodynamic equations from RBE, 
so-called matching conditions corresponding to a desired 
frame are imposed on the higher-order terms from the 
outset, and the hydrodynamic equation in the desired 
frame including Eckart one is formally obtained. It is 
noteworthy that the matching conditions are also a con- 
dition to forbid the appearance of secular terms [10]; we 
refer to Ref. [17] for some problematic aspects of the 
matching conditions. In the present approach, on the 
other hand, we have so far just solved the Boltzmann 
equation in the perturbation theory in a straightforward 
way without imposing any matching conditions, and the 
resulting higher-order terms contain secular terms, which 
apparently invalidates the perturbative expansion for r 
away from the initial time To- 

The key idea in the RG method, however, lies in the 
fact that we can utilize these apparently problematic sec- 
ular terms to obtain an asymptotic solution valid in a 
global domain [19, 20]. Indeed, one can make the follow- 
ing geometrical interpretation of the perturbative solu- 
tion constructed around arbitrary initial time To: That 
is, we have constructed a family of curves / p (t, <t ; To) 
parameterized with To, which curves are all supposed to 
be on the exact solution f p (<r; t) at r = To up to 0(e 3 ), 
although they are admittedly only valid for r near To lo- 
cally. Then, the envelope curve of the family of curves 
will give a global solution in our asymptotic situation, 
which is shown indeed to be the case [19, 20]. Accord- 
ing to the classical theory of envelopes, the envelope that 
is in contact with any curve in the family at t = To is 
obtained [19] by 



_d_ 

dr 



fp(T,a;T ) 



= 0. 



(IV.13) 



The derivative with respect to t hits on the hydro- 
dynamic variables, and hence we have the evolution 
equation of them, which is identified with the hydro- 
dynamic equation [16, 25]. We also note that the in- 
variant manifold [24] corresponding to the hydrodynam- 
ics in the functional space of the distribution function 
is explicitly obtained as an envelope function [16, 17]: 
/e p (t, cr) = f p (r, a; t = t), the explicit form of which 
is referred to Rcf.'s [16, 17]. We note that this solution is 
valid in a global domain of time in the asymptotic region 
[17]. 

Putting back e to 1, Eq.(IV.13) is reduced to the fol- 
lowing form in this approximation, 



(/ P eq + ^ 1) ) = o- 



where 5f p denotes the first-order correction to the dis- 
tribution function 



(i) 



(IV.15) 



If one uses the identity (p ■ u) d/dr + p ■ V = p M <9 M , 
Eq.(IV.14) is found to have the following form 



d„T^ = 0, 



d^Ni* = 0, 



(IV.16) 



withT**" = T^^+ST^ and = N^+SN^. Here, 
the zero-th order terms read T^ ^ = J2 P ^P^P^/p" 3 — 
e/ u v -pA pv and = J2 P jsP^fP = nu p , with e, 

p, and n being the internal energy, pressure, and particle- 
number density for the relativistic ideal gas, respectively, 
while the dissipative parts are given by 



V 



SX P 



(IV. 17) 
(IV.18) 



Note that the dissipative terms are due to the deviation 
dfp 1 ^ of the distribution function from the local one. It is 
well known that the local distribution function only gives 
the (relativistic) Euler equation without dissipation. 



B. Uniqueness of Landau-Lifshitz frame 

In this subsection, we present the dissipative parts 
8T plJ and SN^, explicitly, and discuss their properties. 
An explicit evaluation of Eq.'s (IV. 17) and (IV.18) to- 
gether with (IV.15) gives [16, 17] 



pUff , (IV. 19) 



5T^ = CA^\7-u + 2r 1 A^ f " 7 V p u, 
dX 1 ' = X — V' 1 ^. (IV.20) 



h 2 T 



respectively, with A^ pa = 1/2 ■ (A^A 1 " 7 + A^A^p - 
2/3 • A^A p(y ). Here, h denotes the reduced enthalpy 
per particle. The bulk and shear viscosities and the 
thermal conductivity are denoted by £, r\ and A, respec- 
tively. One readily finds that these formulas completely 
agree with those proposed by Landau and Lifshitz [8]. 
Indeed, the respective dissipative parts ST^ and 5X P 
in Eq.'s (IV.19) and (IV.20) meet Landau and Lifshitz's 
constraints 



5e 
Sn 



= u p ST*" u v = 0, (IV.21) 
= A^6T»»u p = 0, (IV.22) 
= u M SX^ = 0, (IV.23) 



(IV.14) 



which are imposed in a heuristic way in the phenomeno- 
logical derivation [8]. Thus, we have found that the frame 
on which the fluid velocity is defined necessarily becomes 
the Landau-Lifshitz (energy) frame, if the hydrodynam- 
ics is to be consistent with the underlying RBE. 
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Next, let us examine the underlying meaning of Eq.'s 
(IV. 21) to (IV. 23) in terms of the distribution function. 
As was mentioned above, these equations are usually just 
imposed [10] to the higher-order terms of the distribution 
function as the matching conditions without any founda- 
tion to select the hydrodynamic equation in the energy 
frame. We shall clarify that these conditions are equiva- 
lent to the orthogonality condition for the excited modes 
expressed in terms of the inner product [16, 17] and hence 
an inevitable consequence for the relativistic hydrody- 
namics in our analysis which is free from any ansatz. 

We first note that Eq.'s (IV. 17) and (IV. 18) can be 
rewritten as 

ST ^ = J-^p-fp^, (IV. 24) 

p P 

6N " = E/"^p- ( IV - 25 ) 

p 

with 

4> P = (/p eq )- 1 5/i 1) - -[L-'Qoinr'F]^ (IV.26) 

which belongs to the Qo space and thus orthogonal to 
the zero modes, 

(tp%, 4>) =0 for a = 0, 1, 2, 3, 4. (IV.27) 

Recalling the definition Eq. (IV. 10) of the inner prod- 
uct, we see that Eq. (IV.27) with a = /j, = 0, 1, 2, 3 is 
reduced to 

= E^-W^Eiw 

p p 

= u v 5T» v , (IV.28) 

which readily leads to Eq.'s(IV.21) and (IV.22). Quite 
similarly, Eq. (IV.27) with a = 4 is reduced to 

p p 

= u v 5N v , (IV.29) 

which is nothing but Eq.(IV.23). 

We emphasize again that the matching conditions for 
the energy frame are not imposed but uniquely obtained 
in our derivation from the underlying kinetic equation, 
without any assumptions nor ansatz. This facts may 
mean that the relativistic hydrodynamic equation for a 
viscous fluid must be defined in the energy frame, at least 
if it is consistent with the underlying kinetic equation. 

V. SUMMARY AND DISCUSSIONS 

In this paper, we have shown that the renormalization- 
group (RG) derivation of the relativistic dissipative hy- 
drodynamic equation as the infrared dynamics of the un- 
derlying relativistic Boltzmann equation (RBE) uniquely 



leads to the one in the energy frame proposed by Landau 
and Lifshitz, provided that the macroscopic-frame vector, 
which covariantly defines the local rest frame of the fluid 
velocity, is independent of the momenta of constituent 
particles of the system, as it should. 

In relation to other methods of the derivation of hy- 
drodynamic equations based on the RBE, we note that 
we have not assumed any matching conditions [10] but 
uniquely got them in the energy frame from the underly- 
ing kinetic equation, and hence given the foundation for 
the matching conditions. Since any matching conditions 
arc not imposed, secular terms due to the zero modes 
of the linearized collision operator appear inevitably in 
higher-order terms in our approach, but they arc re- 
summed away by the RG/envelope equation to give an 
asymptotic solution valid in a global domain. 

We thus argue that the relativistic hydrodynamic 
equation for viscous fluids must be defined in the energy 
frame, if it is consistent with the underlying kinetic equa- 
tion. Although the RBE which we have adopted as the 
kinetic equation is admittedly suitable only for a dilute 
gas, it is expected that the derived hydrodynamic equa- 
tion itself and hence the uniqueness of the energy frame 
can be valid even for dense systems ; this is found plausi- 
ble if one recalls the universal nature of (non-relativistic) 
Navier-Stokes equation beyond dilute systems, although 
it can be also derived [13, 25] from the (non-relativistic) 
Boltzmann equation. 

Although the present work is confined to the case of the 
so-called first-order equation for a system composed of a 
single component, the uniqueness of the energy frame 
for the relativistic hydrodynamics may keep valid for 
the multi-component systems [10] and the case of the 
so-called second-order causal relativistic hydrodynamics, 
i.e., the extended thermodynamics [26]. In fact, we can 
show that the energy frame is the most natural frame for 
multi-component systems [27] and that the RG derivation 
[28] of the mesoscopic dynamics [30] of the RBE natu- 
rally leads to the extended thermodynamics in the energy 
frame. We can thus assert more firmly that the relativis- 
tic hydrodynamic equations for viscous fluids must be 
defined in the energy frame, if it is consistent with the 
underlying relativistic kinetic equation at all. 
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